Radio-frequency spectroscopy of a strongly interacting spin-orbit coupled Fermi gas 
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We investigate experimentally and theoretically radio-frequency spectroscopy and pairing of a 
spin-orbit-coupled Fermi gas of 40 K atoms near a Feshbach resonance at Bo = 202.2 G. Experi- 
mentally, the integrated spectroscopy is measured, showing characteristic blue and red shifts in the 
atomic and molecular responses, respectively, with increasing spin-orbit coupling. Theoretically, a 
smooth transition from atomic to molecular responses in the momentum-resolved spectroscopy is 
predicted, with a clear signature of anisotropic pairing at and below resonance. Our many-body 
prediction agrees qualitatively well with the observed spectroscopy near the Feshbach resonance. 
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I. INTRODUCTION 

Owing to the unprecedented controllability of inter- 
action and dimensionality, strongly interacting ultracold 
Fermi gases have proven to be an ideal desktop system 
in the study of pairing and superfluidity [l|, [2|. Us- 
ing magnetic field Feshbach resonances, a crossover from 
Bose-Einstein condensates (BECs) to Bardeen-Cooper- 
SchriefTer (BCS) superfluids was successfully demon- 
strated in 2004 |3|, U| and the pairing properties at the 
crossover have been characterized in a number of means 
since then, including particularly radio- frequency (rf) 
spectroscopy [5- 7]. The latest development in this field is 
the realization of a synthetic spin-orbit coupling, which 
couples the pseudo-spin of neutral atoms to their orbital 
motion [§-fl2j|. Such a spin-orbit coupling is responsible 
for a variety of intriguing phenomena in different fields 
of physics. A well-known example is the recently discov- 
ered topological insulators in solid-state pjl 03] • In the 
context of ultracold atomic Fermi gases, it is therefore 
natural to ask: what is the consequence of the interplay 
of strong interaction and spin-orbit coupling? 

In this paper, we investigate rf-spectroscopy of a 
strongly interacting spin-orbit coupled Fermi gas of 40 K 
atoms. The previous works on spin-orbit coupled Fermi 
gas explored essentially the non-interacting limit flll[l2l|. 
The current work is the first demonstration of effects of 
spin-orbit coupling in an interacting Fermi gas. In recent 
BEC-BCS experiments, rf-spectroscopy has been particu- 
larly useful in studying pairing and superfluidity, yielding 
information about the pairing gap [5j and pair size [6j. 
Furthermore, momentum-resolved rf-spectroscopy gives a 
direct information of the low-energy excitation spectrum 
and quasiparticles [7]. Here, by developing a many-body 
T-matrix theory we show theoretically that both atomic 
and molecular responses in the rf-spectroscopy, aris- 
ing respectively from free fermionic atoms and bosonic 



molecules, are greatly modified by spin-orbit coupling. 
In particular, the resulting anisotropic pairing, domi- 
nated by the two-body effect on the BEC side of the 
Feshbach resonance and by the many-body effect near 
resonance, is clearly evident in the momentum-resolved 
spectroscopy. Experimentally, we measure the integrated 
rf-spectroscopy and report characteristic blue and red 
shifts in the atomic and molecular responses, respectively, 
which are in good agreement with theory. 

The remainder of this paper is organized as follows. 
In Sec. II, we describe briefly the experimental setup 
and the model Hamiltonian. In Sec. Ill, we present the 
experimental and theoretical results of radio-frequency 
spectroscopy near Feshbach resonances. We introduce 
briefly a many-body T-matrix theory in Sec. III(A), and 
in Sec. III(B) and Sec. III(C) we discuss respectively 
the integrated spectroscopy of bound molecules and the 
momentum-resolved spectroscopy of fermionic pairs at 
resonance. The comparison between experiment and the- 
ory is reported in Sec. III(D). Finally, we conclude in Sec. 
IV. The Appendix A is devoted to solving the many-body 
T-matrix theory within the pseudogap approximation. 



II. EXPERIMENTAL SETUP AND MODEL 
HAMILTONIAN 

The experimental setup has been described in our pre- 
vious works fill , [l5l - [l7| , in which a Bose- Fermi mixture 
of 40 K and 87 Rb atoms is cooled using well-developed 
evaporative and sympathetic cooling techniques in a 
quadrupole-Ioffe configuration magnetic trap, and is 
transported into an optical trap. A degenerate Fermi gas 
of about N ~ 2 x 10 6 40 K atoms in the \F = 9/2, m F = 
9/2) internal state is evaporatively cooled to tempera- 
ture T/Tp ~ 0.3 with bosonic 87 Rb atoms, where Tp 
is the Fermi temperature defined by Tp = Ep/ks = 
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FIG. 1: (color online) (a) and (b) Experimental realization 
of a strongly-interacting Fermi gas of 40 K atoms with spin- 
orbit coupling, (c) The integrated rf-spectroscopy below the 
Feshbach resonance (at B = 201.6 G and a s ~ 2215. Gas, 
where as is the Bohr radius), in the presence (solid cir- 
cles) and absence (empty circles) of the spin-orbit coupling. 
The Raman detuning is 5 — 0. The dimensionless interac- 
tion parameter 1/(/cfci s ) ~ 0.66. The fraction is defined as 
^-5/2/(^-5/2+^-7/2), where N_ 5/2 and N_ 7/2 are obtained 
from the TOF absorption image. SOC: spin-orbit coupling. 



(QNy^hO/kB and uJ c± 2tt x 130 Hz is the geometric 
mean trapping frequency. 87 Rb atoms in the mixture are 
then removed by a 780 nm laser pulse. Subsequently, 
fermionic atoms are transferred into the lowest hyperfine 
state \F = 9/2, mi? = —9/2) via a multi-photon rapid 
adiabatic passage induced by a radio frequency field at 
lower magnetic field. To prepare a two-component 40 K 
Fermi gas in an equal mixture of |t) = \F = 9/2, mp = 
-7/2) and ||) = \F = 9/2, m F = -9/2) states, a homo- 
geneous bias magnetic field, produced by the quadrupole 
coils (operating in the Helmholtz configuration), is raised 
to about B « 219.4 G and then a radio frequency ramp 
around 47.45 MHz is applied for 50 ms. To create strong 
interactions, the bias field is ramped from 204 G to a 
value near the Bo = 202.2 G Feshbach resonance at a 
rate of about 0.08 G/ms. 

We create spin-orbit coupling using Raman process 
[slGj]]? as illustrated in Figs. 1(a) and 1(b). A pair of Ra- 
man beams from a Ti-sapphire laser counter-propagate 
along the x axis and couple the two spin states. The 
intensity of beams is I = 50 mW and their frequencies 
are shifted respectively by 75 and 120 MHz, using two 
single-pass acousto-optic modulators. These two Raman 



beams intersect in the atomic cloud with 1/e 2 radii of 
200 jam and are linearly polarized along z and y axis 
directions, respectively. The momentum transferred to 
atoms during the Raman process is |q^| = 2k r sin(0/2), 
where k R = 2tt/\r is the single-photon recoil momen- 
tum, Xr = 772.1 nm is the wavelength, and = 180° is 
the intersecting angle of two Raman beams. In the sec- 
ond quantization, this Raman process may be described 
by the term 



H R = 



J dr (r) e 1 



2k R x 



<Mr)+H.c. , (1) 



where (r) is the creation field operator for atoms in 
the spin-state a =t, | and Qr is the coupling strength of 
Raman beams. For a detailed discussion on the Raman 
coupling strength Qr, we refer to the recent theoretical 
work by Wei and Mueller [18] . 

In this paper, we use a larger bias magnetic field than 
the one used in our previous study [11], in order to cre- 
ate strong interactions. Due to a decoupling of the nu- 
clear and electronic spins, the Raman coupling strength 
decreases with increasing the bias field [18]. To compen- 
sate this reduction, here we use a smaller detuning of the 
Raman beams with respect to the atomic "Dl" transition. 

To see clearly the spin-orbit coupling in our setup, it 
is convenient to take a gauge transformation, ^t( r ) = 
e lkRX ^^(r) and V^( r ) = e~ %kRX ^ \(v). Our system may 
therefore be described by a model Hamiltonian H = Ho + 
Hint, where 



Ho 



H 2 (k±k R e x ) 
r) ^ —^*(r) + 



dr 



2m 

tf J (r) ^ (r) + H.c 



is the single-particle Hamiltonian and 

-Hint = U J dr*| (r) (r) (r) * t (r) 



(2) 



(3) 



describes the contact interaction. 
Eq. ([2]) we have used k = zV, u - 



In the first line of 
" for a =t and "-" 



for a Using the Pauli matrices a x: a y and <r 2 , the 
single-particle Hamiltonian may be rewritten in a com- 
pact form, 



Ho 



I 



dr& 



2m 



hcr x + \k x a z 



<*>, (4) 



where the spinor field operator <E>(r) = [\E^ (r) , (r)] T . 
We have defined a spin-orbit coupling constant A = 
h 2 kR/m and an "effective" Zeeman field h = Qr/2. 

To create a strongly interacting Fermi gas with spin- 
orbit coupling, after the bias magnetic field is tuned to a 
final value B (which is varied) , we ramp up adiabatically 
the Raman coupling strength in 15 ms from zero to its fi- 
nal value = 1.5 Er with Raman detuning 5 = 0, where 



3 



the recoil energy E R = H 2 k 2 R /(2m) ~hx 8.36 kHz. The 
temperature of the Fermi cloud after switching on the 
Raman beams is at about 0.67>, as in our previous mea- 
surement for a non-interacting spin-orbit coupled Fermi 
gas [jjj. The Fermi energy are Ep — 2.§Er and the 
corresponding Fermi wave vector is hp — 



III. RADIO FREQUENCY SPECTROSCOPY 

To characterize the strongly interacting spin-orbit cou- 
pled Fermi system, we apply a Gaussian shape pulse of rf 
field with a duration time about 400 /is and frequency uo 
to transfer the spin-up fermions to an un-occupied third 
hyperfine state |3) = \F = 9/2, m F = -5/2). The Gaus- 
sian shape pulse is generated by the voltage-controlled rf 
attenuators. The Gaussian envelope hence results in the 
elimination of the side lobes in rf spectra. The Hamilto- 
nian for rf-coupling may be written as, 



Vrf = V dr 



4 (r) * t (r) + H.c. 



(5) 



where (r) is the field operator which creates an atom 
in 1 3) and Vb is the strength of the rf drive. The effective 
momentum transfer kue x in V r f results from the gauge 
transformation. After the rf pulse, we abruptly turn off 
the optical trap, the magnetic field and the Raman laser 
beams, and let the atoms ballistically expand for 12 ms in 
a magnetic field gradient applied along z and take time- 
of-flight (TOF) absorption image along y. We measure 
the spin population of the final state |3) for different rf 
frequencies to obtain the rf spectra T(uj). 

For a weak rf drive, the number of transferred fermions 
can be calculated using linear response theory. At this 
point, it is important to note that the final-state interac- 
tions for 40 K atoms in the third state and in the spin-up 
or spin-down state is typically small Theoretically, 
the rf transfer strength at a given momentum is there- 
fore determined entirely by the single-particle spectral 
function of spin- up atoms A^: 



T(k,cj) = Aft(k + k R e x ,6k - ii - 
/(ek - M - + &*>3t)j 



(6) 



where ek = h 2 k 2 / (2m), \i is the chemical potential of the 
spin-orbit system, hw^ c± h x 47.1 MHz is the energy 
splitting between the third state and the spin-up state, 
f(x) = \j(e x l kBT + 1) is the Fermi distribution function, 
and we have taken the coupling strength Vo = 1. For 6 Li 
atoms, however, the final state effect is usually significant 
[6(| . The rf transfer strength will no longer be simply 
determined by the single-particle spectral function. 

Experimentally, one could measure the momentum- 
resolved rf spectroscopy along the ^-direction T(k x ,uj) = 
J2k y k z r(k,cj), or, after integration obtain the fully inte- 
grated rf spectroscopy T(uj) = T(k, uo). Due to small 
signal to noise ratio, we currently have difficulty obtain- 
ing momentum-resolved rf signal experimentally. 



In Fig. 1(c), we show the integrated rf-spectroscopy of 
an interacting Fermi gas below the Feshbach resonance, 
with or without spin-orbit coupling. Here, we carefully 
choose the one photon detuning of the Raman laser to 
avoid shifting Feshbach resonance by the Raman laser 
on the bound-to-bound transition between the ground 
Feshbach molecular state and the electronically excited 
molecular state. The narrow and broad peaks in the 
spectroscopy should be interpreted respectively as the 
rf-response from free atoms and fermionic pairs. With 
spin-orbit coupling, we find a systematic blue shift in the 
atomic response and a red shift in the pair response. The 
latter is an unambiguous indication that the properties 
of fermionic pairs are strongly affected by spin-orbit cou- 
pling. 



A. Many-body T-matrix theory 

Let us now consider theoretical understanding of the 
observed red shift for fermionic pairs. Near Feshbach 
resonances, it is important to treat atoms and fermionic 
pairs on an equal footing. For this purpose, it is con- 
venient to develop a many-body theory within the T- 
matrix approximation by summing all ladder diagrams 
fl9l [20j. In the presence of spin-orbit coupling, it is nec- 
essary to define a finite-temperature Green function 



g(r ) r';r>0) = -($(r,r) $t ( r ', )> , 



(7) 



which is a 2 by 2 matrix even in the normal state. We 
adopt a partially self-consistent T-matrix scheme and 
take one non-interacting and one fully dressed Green 
function in the ladder diagrams [20|. The summation 
of all ladder diagrams leads to the Dyson equation, 



Q{K)=[Qz\K)-V{K)]\ 
where the self-energy is given by 

E(K) = J2 [t(Q)(iv y )Go(K - Q)(ia y ) 



Here 



t(Q) 



Uo 



(8) 



(9) 



(10) 



i + Uox(Q) 

is the (scalar) T-matrix with a two-particle propagator 



X 



(Q) = 1 £ TV \g(K) (ia y ) G (K - Q) (ia v )} (11) 

1 K L J 



and the non-interacting Green function 
Go(K) = [iu) m - e k + M - Er ~ ha x - Xk x a z 



(12) 



For convenience, we have used the short notations K 
(k,zu; m ), Q = (q,i^ n ), and J2k = ^ T X)ku; > and u 
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FIG. 2: (color online) Evolution of the predicted integrated rf- 
spectroscopy as a function of the Raman coupling strength. 
Here, l/k F a s = 0.66, T = 0.6T F , and k F = 1.6k R . The 
inset shows a linear contour plot of the predcited momentum- 
resolved spectroscopy at Qr — 2Er. 




and v n are respectively the fermionic and bosonic Mat- 
subara frequencies. We have also defined a Green func- 
tion for holes, Q(K) = - [G(-K)} T . 

Eqs. (|8l)- ([T2|) generalize the earlier T-matrix diagram- 
matic theory without spin-orbit coupling [20] . Here the 
modification arising from the spin-orbit coupling includes 
the use of Green functions in a general 2 by 2 matrix form 
and accordingly the appearance of the vertex ia y in Eqs. 
(O and (fTTj) . The derivation of these T-matrix equa- 
tions is too technical and is not the focus of this work. 
Therefore, we will discuss the details of the derivation 
elsewhere. 

In general, the self-consistent T-matrix equations are 
numerically difficult to solve. At a qualitative level, how- 
ever, we may adopt a pseudogap decomposition for the 
T-matrix [2l[ and obtain a set of coupled equations for 
the chemical potential and pairing gap. In the Appendix 
A, we discuss in detail the pseudogap approximation. By 
solving the coupled equations (see, i.e., Eqs. (|A8[) and 
(|A9P ), we calculate the single-particle spectral function 



(13) 



and the rf-transition strengths T(k x1 uo) and T(uo). To 
take into account the experimental energy resolution of 
the spectroscopy 7 ~ 0.1 Er [12|, we replace the Dirac 
delta function (which appears in *4^(k, uj)) by 



&{x) 



7/tt 
x 2 + 7 



2 ' 



(14) 



B. Integrated spectroscopy of bound molecules 

In Fig. 2, we predict the integrated rf-spectroscopy 
of a spin-orbit coupled Fermi gas at B = 201.6 G and 



FIG. 3: (color online) (a) Predicted phase diagram of a 
strongly-interacting spin-orbit coupled Fermi gas at Qr = 
1.5Er and k F — (b)-(d) Linear contour plots of the 

zero temperature momentum-resolved rf-spectroscopy across 
the Feshbach resonance, in arbitrary units. The inset in 
(a) shows the momentum-resolved rf spectroscopy at the 
pseudgap temperature T* in the unitary limit. In theoret- 
ical calculations, we have set uj^ — 0. 



T = 0.6Tp. With increasing the strength of Raman 
beams from to 4Er, the atomic and pair peaks shift to 
the higher and lower frequencies, respectively, in qualita- 
tive agreement with the experimental observation in Fig. 
1(c). At this interaction parameter (l/(/ci?a s ) ~ 0.66), 
the red shift of pair peaks may be understood from 
the binding energy of pairs in the two-body limit: the 
stronger effective Zeeman field ft, the smaller binding 
energy of two-particle bound states [22] . In the inset, 
we show the prediction of the momentum-resolved rf- 
spectroscopy at = 2Er. It is highly asymmetric as a 
function of momentum. The predicted atomic response is 
in good agreement with the experimental observation for 
a non- interacting spin-orbit-coupled Fermi gas fill [l2|. 



Momentum resolved spectrosocpy near 
Feshbach resonances 



We now turn to the rf-spectroscopy in the vicinity of 
the Feshbach resonance. In Fig. 3(a), we plot the super- 
fluid transition temperature T c and the pairing break- 
ing (pseudogap) temperature T* of a spin-orbit coupled 
Fermi gas at Qr = 1.5i^R and kp = 1.6/c#. The pseu- 
dogap temperature is calculated using the standard BCS 
mean-field theory without taking into account the pre- 
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formed pairs (i.e., A pg = 0) [21]. We find that the re- 
gion of superfluid phase is suppressed by spin-orbit cou- 
pling. In particular, at resonance the superfluid transi- 
tion temperature is about T c ~ 0. 1297V, smaller than 
the measured value of T c ~ 0.167(13)Tp [23] or the pre- 
dicted value of T c ~ 0.1 5Tp (under the same pseudogap 
approximation) for a unitary Fermi gas in the absence 
of spin-orbit coupling. Thus, experimentally it becomes 
more challenge to observe a spin-orbit coupled fermionic 
superfluid in the present experimental scheme. 

In Figs. 3(b)-3(d), we show the zero-temperature 
momentum-resolved rf-spectroscopy across the reso- 
nance. On the BCS side (l/fc^c^ — — 0.5), the spec- 
troscopy is dominated by the response from atoms and 
shows a characteristic high-frequency tail at k x < 
[H d [H|. Towards the BEC limit (l/k F a s = 0.5), 
the spectroscopy may be understood from the picture of 
well-defined pairs and shows a clear two-fold anisotropic 
distribution [25|. The spectroscopy at the resonance 
is complicated and should be attributed to many-body 
fermionic pairs. The change of spectroscopy across reso- 
nance is continuous, in accord with a smooth BEC-BCS 
crossover [2|]. 

In the inset of Fig. 3(a), we show the momentum- 
resolved rf-spectroscopy at the resonance and at the pseu- 
dogap pairing temperature T*. It is interesting that the 
anisotropic distribution survives well above the super- 
fluid transition temperature T c . An experimental obser- 
vation of such a spectroscopy would be a strong indica- 
tion of the anisotropic pseudogap pairing of a spin-orbit 
coupled Fermi gas in its normal state. 

D. Experiment vs Theory 

In Fig. 4, we examine our many-body theory by com- 
paring the theoretical predictions with the experimental 
data for integrated rf-spectroscopy. At the qualitative 
level, we do not consider the trap effect and take the rel- 
evant experimental parameters at the trap center. Oth- 
erwise, there are no adjustable free parameters used in 
the theoretical calculations. As shown in Fig. 4, we find 
a qualitative agreement between theory and experiment, 
both of which show the red shift of the response from 
fermionic pairs, as we already discuss in Sec. III(B). The 
worst agreement occurs close to the Feshbach resonance 
(Fig. 4(c)), where our many-body theory fails to take into 
account properly the strong interactions between atoms 
and pairs. 



IV. CONCLUSION 

In conclusion, we have investigated experimentally and 
theoretically rf-spectroscopy and fermionic pairing in a 
strongly-interacting spin-orbit coupled Fermi gas of 40 K 
atoms near a Feshbach resonance. A red shift of the re- 
sponse from fermionic pairs, induced by spin-orbit cou- 



pling, is observed in integrated rf-spectroscopy below the 
resonance, in qualitative agreement with a many-body T- 
matrix calculation. Momentum-resolved rf-spectroscopy 
of fermionic pairs has been predicted across the resonance 
at all temperatures, showing a characteristic anisotropic 
distribution. This is to be confronted in future experi- 
ments. 

We note that the typical experimental temperature in 
this work is about 0.61>. In the future, we wish to re- 
duce the temperature of the strongly interacting Fermi 
cloud of 40 K atoms down to 0.2T^, close to the superfluid 
transition. This is encourged by a recent calculation by 
Wei and Mueller [18], who showed that the heating of 
the cloud due to Raman transition is not significantly 
affected by the magnetic field needed for Feshbach reso- 
nances. 

Theoretically, it is possible to solve the many-body 
T-matrix theory without the pseudogap approximation. 
On the other hand, at the relatively high temperature 
T ~ 0.6Tp, alternatively we may use a virial expansion 
theory to obtain quantitative predictions for the radio- 
frequency spectroscopy in harmonic traps [271-129] . These 
possibilities will be addressed in later studies. 
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Appendix A: Pseudogap approximation 

The pseudogap approximation is advanced by the 
Chicago group [21]. In this approximation the T-matrix 
is separated into two parts, t(Q) = t sc (Q) + t pg (Q), so 
that the contribution from the superfluid order parame- 
ter for condensed pairs A 2 ^, 

t sc (Q) = -^5(Q), (Al) 

and the contribution from the pseudogap for un- 
condensed pairs, 

^ = -Eu«). ( A2 ) 

become explicit. The full pairing order parameter is given 
by A 2 = A 2 C + A 2 . Accordingly, we have the self-energy 
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FIG. 4: (color online) Comparison between theory and experiment for the integrated rf-spectroscopy. In theoretical simulations, 
we use Qr = 1.5Er, Jcf = and T = 0.6Tf, according to the experimental setup. The solid circles (lines) and empty 

circles (dashed lines) show respectively the experimental data (theoretical predictions) in the presence and absence of spin-orbit 
coupling with Raman detuning 6 = 0. The dimensionless interaction parameter l/(fea s ) in (a), (b), and (c) are 0.89, 0.66, 
and 0.32, respectively. 



H(K) = ^sc(K) + ^p g (K) (IH, where 

S sc = -A 2 sc (ia y )Go(K)(ia y ) 

and 



(A3) 
(A4) 



To obtain Yi pg in the above equation, it is assumed that 
the pair propagator x(Q) peaks around Q = [2l|. 



We note that, at zero temperature the pseudogap ap- 
proximation is simply the standard mean- field BCS the- 
ory, in which 



X(K) = -A 2 (ia y )g (K)(ia y ). 



(A5) 



Above the superfluid transition, however, it captures the 
essential physics of fermionic pairing and therefore should 
be regarded as an improved theory beyond mean-field. 
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To calculate the pseudogap A pg , we approximate 

tpg(Q -0) = Z [iv n - ^ q + flpair] , (A6) 

where the residue Z and the effective dispersion of pairs 
^ q = h 2 q 2 /2M* are to be determined by expanding 

X {Q) about Q = [Ulii]. The form of tpgiQ) leads 
to 

A 2 pg (T) = Z- 1 /B(«q - Impair), (A7) 

q 

where = \j{e x l kBT — 1) is the bosonic distribution 

function. We arrive finally at two coupled equations, the 



gap equation 

7T +X(<2 = 0) = Z^r (A8) 
and the number equation 

n = ^TrG(K), (A9) 

K 

to determine the superfluid order parameter A sc and the 
chemical potential /i, respectively, for a given set of pa- 
rameters. 
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